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A novel point of view on the phenomenon of self-pulsations is presented, which shows that they are 
a balanced state formed by two counteracting processes: beating of modes and bistable switching. 
A structure based on two coupled nonlinear microcavities provides a generic example of system 
with enhanced abihty to this phenomenon. The specific design of such structure in the form of 
multilayered media is proposed, and the coupled mode theory is applied to describe its dynamical 
properties. It is emphasized that the frequency of self-pulsations is related to the frequency splitting 
between resonant modes and can be adjusted over a broad range. 
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I. INTRODUCTION 

Optical bistability is a general nonlinear phenomenon, 
which can be observed in a large variety of systems pro- 
vided that nonlinearity is complemented by some form 
of feedback. Apart from the classical Fabry-Perot res- 
onators, the possible configurations of bistable systems 
include multilayered structures [2 [5] , fiber Bragg grat- 
ings |31 13] and microcavities in photonic crystal waveg- 
uides 016]. 

It was noticed however that the description of these 
systems in the frequency domain does not capture the full 
complexity of their temporal behavior. After overcoming 
the first switching threshold, the transient response of 
these systems may change dramatically, leading to self- 
pulsations or even chaos [7J{5j. These instabilities in the 
time domain are usually explained from the viewpoint of 
nonlinear dynamics as a manifestation of the Hopf bifur- 
cations [TUl414j . Although such description is appropri- 
ate, it is often desirable to have a control over these insta- 
bilities, as they might be useful for certain applications. 
In this work, we present a novel point of view on the 
phenomenon of self-pulsations which shows that they are 
a balanced state formed by two counteracting processes: 
beating of modes and bistable switching. We apply this 
concept to design a structure where self-pulsations can 
be excited very efficiently. 

The paper is organized as follows. In Section [nj we 
consider a typical example of the structure that supports 
self-pulsations. It is represented by the Bragg grating, 
and we explain how the beating of modes and bistable 
switching can interplay in this particular case. In Sec- 
tion |IIH we propose a design of multilayered structure. 



which is based on two coupled nonlinear microcavities 
and shows pronounced self-pulsations. The structure has 
several independent parameters that make the adjust- 
ment of its properties easier. In Section |IV[ we derive a 
set of coupled mode equations which describes the switch- 
ing dynamics in such system and discuss how the qual- 
itative behavior of the structure agrees with the results 
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obtained by the linear stability analysis. As a possible 
application, we show that self-pulsations can be used to 
convert a continuous wave signal into a regular train of 
ultrashort pulses. The last section summarizes the re- 
sults and presents the conclusions. 



II. BEATING OF MODES 

As was mentioned in Introduction, self-pulsations tend 
to appear after the first switching threshold on the hys- 
teresis curve. However, each loop of the hysteresis curve 
is related to a resonance in the linear transmission spec- 
trum. For example, a Bragg structure obtained by pe- 
riodic alternation of two different layers is considered 
on Fig. [l] It has many resonances near the band edge 
[Fig. [T]^a)] which can be bent into the band gap region 
by the Kerr nonlinearity producing a multistable hys- 
teresis curve [Fig. [ijb)] . Each resonance creates its own 
loop when input intensity is gradually increased and de- 
creased. The difference between the first and higher 
switching thresholds is clearly pronounced in the field 
profiles of corresponding resonances. It turns out that 
the first resonance (counting from the upper band edge) 
has only a single localization center in the field profile, 
the second has two such centers [Fig. [l|c)] , and this num- 
ber continues to increase by one for the next resonances 
[TSHTT] . Therefore, even the Bragg structure without any 
defects can be considered as a set of coupled microcavi- 
ties. 

It is known however that the beating of modes can 
exist in coupled systems. Although it is a purely linear 
effect, it can so strongly redistribute the energy between 
constituent parts of the system that they will meet condi- 
tions for up- and down-switching periodically. The need 
to overcome the first switching threshold in the case of 
Bragg structures is thus related to the fact that they 
can be approximated as a single-mode microcavity near 
the first resonance, and the beating of modes is not pos- 
sible in such circumstances. When the input intensity 
reaches the switching threshold, the higher order reso- 
nances shifted by the Kerr nonlinearity will dominate, 
and the beating of modes can be excited. In principle. 
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FIG. 1. (Color online), (a) Linear transmission spectrum 
of the Bragg structure consisting of 64 layers. For conve- 
nience, the resonances are counted from the upper band edge 
and marked by the red triangles with numbers. The posi- 
tive Kerr nonlinearity can bend the linear resonances to the 
lower frequencies, (b) A typical hysteresis curve for signals 
with carrier frequency taken inside the band gap region. The 
loops created by the resonances for the increasing and de- 
creasing input intensity are shown by arrows, (c) The linear 
field profile of the second resonance with one of the localiza- 
tion centers highlighted by the red rectangle. The background 
shows alternation of layers inside the structure. 



self-pulsations can be observed even without overcoming 
any threshold, if one chooses the carrier frequency to be 
initially near a higher order resonance. 

There are two requirements for the onset of beating. 
Firstly, the system should support normal modes with 
slightly different frequencies. Most often, they can be 
built as symmetric and antisymmetric combinations of 
a degenerate mode. The difference of their frequencies 
determines the period of the beating. Secondly, the sys- 
tem should be excited in a way that creates an initial 
imbalance of energy between its parts, for example, by 
rapidly varying input intensity. The latter requirement 
is particularly important for open systems, because beat- 
ing as any other transitional response is prone to decay 
in time. If the initial imbalance of energy is not regu- 
larly replenished, the system will settle in a stationary 
state. It turns out however that the systems of coupled 
nonlinear microcavities are able to sustain the beating of 
modes for arbitrary long time even if the input intensity 
is constant. This happens naturally because switching it- 



self changes energy stored in microcavities in a step-like 
manner. 

Therefore, switching and beating are closely related 
processes. Each of them can excite the other, and at 
the same time it needs to be sustained by the other. 
The interplay between these processes can take a peri- 
odic form, resulting in an infinite series of switchings, or 
self-pulsations. 

III. COUPLED MICROCAVITIES 

The analysis performed in the previous section leads 
to conclusion that the simplest system with the ability 
to self-pulsations should consist of just two microcavities, 
and each of them can be strictly bistable. We consider 
a possible design of such system in the form of multilay- 
ered structure. It is sufficient to have dielectric layers of 
two types, which will be denoted as 'L' and 'H'. They 
correspond to materials with different refractive indices, 
and it will be assumed that the optical thickness of these 
layers satisfies the quarter wave condition. In what fol- 
lows, the specific materials used for the layers 'L' and 
'H' are polydiacetylene 9-BCMU with linear (nonlinear) 
refractive index riL — 1.55 (n2L — 2.5 • lO^'^ cm^/MW) 
and rutile with nn = 2.3 (nan = 10"^ cmVMW), re- 
spectively [THJ [H] . The quarter wave condition is set to 
Ao = 0.7 ^m, which gives the thicknesses of the layers 
c^L — 112 nm and dn = 76 nm. 

The microcavities can be formed by breaking period- 
icity in the arrangement of layers. We propose the fol- 
lowing symbolic formula for the multilayered structure: 

(HL)P+™(LH)"L(HL)"(LH)™+'', (1) 

where the parameters m and p are integer numbers which 
satisfy the conditions m > 1 and p> 0. It is important to 
consider their meaning in more detail and to emphasize 
how they are related to various functional parts of the 
structure [Fig. [2]ja)]. 

The two major parts are represented by the combina- 
tion (IIL)™(LII)™, which can be viewed as a result of im- 
perfect junction of two periodic sequences. This causes 
the formation of a point-like defect in the middle, and 
creates a strong resonance in the transmission spectrum. 
This resonance is located exactly at the frequency for 
which the quarter wave condition holds, and the mirror 
symmetry of the structure ensures that it demonstrates 
perfect transmission [Fig. ^h)]. The quality factor Q 
of this resonance can be adjusted over a broad range by 
varying the number m. When the refractive indices of 
the layers riL and nn have a large contrast, the following 
formula can be used to estimate the quality factor 

Q-^-jP^J'^r. (2) 

27 4(nH - riL) \nL J 

The purpose of the single layer 'L' placed in the middle 
of the structure is twofold. Firstly, it smoothly connects 
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FIG. 2. (Color online), (a) Example of multilayered structure 
obtained from the symbolic formula with various func- 
tional parts highlighted, (b) Transmission spectrum of the 
single microcavity for m = 3 (shown as an inset on the left) 
that was computed with the transfer matrix method (solid 
blue) and the coupled mode theory (dashed red). The res- 
onance of perfect transmission caused by the point defect is 
in the middle of the band gap region, where the two spectra 
overlap, (c) Transmission spectrum of two coupled microcavi- 
ties for m = 3 and p = 0. Due to the coupling, the degenerate 
resonance splits into a doublet of even and odd modes. 



the two microcavities without introducing an additional 
defect into the structure. Secondly, it changes the phase 
mismatch between microcavities. If it were absent, the 
interaction between them would be suppressed, and the 
excitation of beating would not be possible. In addition, 
the structure can be extended with Bragg mirrors of the 
form (IIL)P and (LH)'' on the sides to make internal pro- 
cesses less affected by too early leakage of energy. 

The influence of the coupling manifests itself already 
in the linear transmission spectrum, where the otherwise 
degenerate resonances of two identical microcavities split 
into a doublet [Fig.[2|^c)]. Similar to the bonding and an- 
tibonding states of diatomic molecules, it is possible to 
distinguish the symmetric (w < loq) and antisymmetric 
(w > Wo) modes in this doublet and to expect the onset of 
beating, when energy periodically tunnels from one mi- 
crocavity to another. If the structure is nonlinear, these 
oscillations can modulate the optical properties leading 
to Raman-like scattering of incident waves on it and sclf- 
pulsations. 



IV. SWITCHING AND SELF-PULSING 

To describe the switching dynamics in a single micro- 
cavity, we apply the coupled mode theory (CMT), which 



has been recognized as a useful tool for studying the prop- 
erties of linear and nonlinear defects in photonic crystal 
waveguides [5DH22]. Assuming that the microcavity is 
mirror symmetric, the CMT equations [see Fig. |3]^a) for 
notation] can be written as 
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where A{t) is the amplitude of the mode Eq{x) with res- 
onant frequency luq and damping constant 7, u±{v±) are 
the amplitudes of in- and out-going plane waves on the 
left (right) side of the cavity. 

Although the materials are assumed to be lossless, 
damping takes place due to the leakage of energy out 
of the structure and the decay rate can be estimated as 
7 = [(1/c) J^^ e\EQ\'^dx]~^. The integration in this for- 
mula is performed over the full length of the microcavity 
(spanning from to L) , c is the speed of light in the vac- 
uum, and e{x) is the linear permittivity (e — "n?). Since 
the electric field inside the microcavity is a product of 
the mode amplitude A{t) and the mode profile Eo{x)^ it 
is convenient to consider the mode profile as a dimen- 
sionless function which satisfies the boundary conditions 
Eq{0) = Eq{L) = 1. This ensures that the expression for 
the damping constant has a proper dimensionality. 

The nonlinearity is responsible for the shift of reso- 
nant frequencies and can be taken into account by in- 
troducing the characteristic intensity for the microcavity 
/q = [(loq/c) n2e\EQ\'^dx]~^ . This formula follows from 
the orthogonal projection of the perturbation term which 
contains the nonlinear index of refraction 712(2;) on the 
microcavity mode. 

For the structure described by m = 3, the damping 
constant and characteristic intensity take the following 
values 7"^ = 27.6 fs, Iq = 154 MW/cm^. The question 
how accurate the CMT model is can be answered by com- 
paring the linear transmission spectrum computed with 
a rigorous transfer matrix method and with the following 
formula derived from Eqs. ([3])-(|4| 
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where /in = (/out = I'l'-P) is the intensity of 

incident (transmitted) wave. Despite of the fact the 
quality factor of the considered resonance is rather low 
(Q = 37.2), the two spectra overlap in a wide range of 
frequencies around ujq [Fig. [2jb)]. 

For a fixed value of frequency detuning, Eq. ([s]) de- 
scribes the hysteresis curve of the single microcavity 
[Fig. |3jb)]. It is strictly bistable and the switching re- 
gion can be determined from the linear stability analy- 
sis. The points on hysteresis curve, which have at least 
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FIG. 3. (Color online), (a) A sketch of the model used to de- 
scribe coupling between microcavities. The response of each 
element is approximated by a scattering matrix, (b) The hys- 
teresis of output intensity for a single microcavity. The real 
and imaginary parts of eigenfrequencies obtained from the 
linear stability analysis are superimposed on the right side. 
Switching occurs when Im(a;p) > and Re(ajp) = 0. (c) The 
hysteresis of output intensity for two coupled microcavities. 
Self-pulsing occurs when Im(a;p) > and Ke{ujp) ^ 0. 



one eigenfrequency with positive imaginary part, cannot 
be stable, because small perturbations will exponentially 
grow forcing switching to another branch of hysteresis. 
It is worth noting that the eigenfrequencies of the lin- 
earized system ujp should satisfy a quadratic equation, 
{ujp/j+i)'^+q = 0, where g is a real coefficient. Therefore, 
the transient response after the switching should have the 
form of exponentially decaying oscillations exp(— iwpi), 
where Wp/7 = iby^lqj — i. 

A qualitatively different behavior is possible in pres- 
ence of two coupled microcavities. The corresponding 
CMT equations can be written as 
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where til is the refractive index of the phase-matching 
layer 'L', and the parameter k — {—nii/ni^)P takes into 
account the influence of the Bragg mirrors. The linear 
normal modes of the system ([6])-([7]) are given by com- 
binations Ai ± A2, and their resonant frequencies are 
w = Wo T ni^l/'^- It is important that the frequency 



FIG. 4. (Golor online), (a) The intensities of incident and 
transmitted signals as functions of time during the switch- 
ing between stable states and after transition to the self- 
pulsations, (b) A phase portrait of the switching, which is ob- 
tained from the real and imaginary parts of transmission am- 
plitude. The final state corresponds to a stable spiral point, 
(c) A similar phase portrait for one period of self-pulsations. 
A limit cycle is formed due to the Hopf bifurcation [12j . 



splitting depends only on the number m and can be es- 
timated using Eq. ([2]), while the number p determines 
the external quality factor Qcxt — k^Q and can lower the 
effective characteristic intensity IcS — lo/^'^- The latter 
property can be useful if materials have a weak Kerr non- 
linearity, but for simplicity it will be assumed that p = 
and consequently k = 1. 

The hysteresis curve of this system shows that there is 
a range of input intensities where only unstable solutions 
are possible [vertical gray region in Fig.jsj^c)]. This range 
falls exactly between the resonances of the doublet, and 
the transient dynamics in this region is different from the 
switching [see Fig. |4] and caption for details]. In the time 
domain simulation, the input intensity of a continuous 
wave signal was slowly varied to examine stability of sev- 
eral key points on the hysteresis curve. After overcoming 
the second threshold the system goes into a state of in- 
finite switching. The period of these pulses extracted 
from the simulation data was 252 fs, which is of the same 
order of magnitude as the beating period of the two reso- 
nant modes 47r/(nL7) = 224 fs. Therefore, the frequency 
of self-pulsations cogp — 4 THz can be estimated by the 
following formula 



Wsp = (nLWo)/(4Q). 



(8) 



It shows explicitly that smaller internal quality factors 
can produce train of pulses with higher repetition rate. 
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since photons tend to spend less time in each cavity and 
they exchange energy at a faster rate. This property can 
be useful in optical communications, where such pulses 
can transfer bits of information [55^. Another possible 
application is optical clocks for photonic circuits, where 
train of pulses can synchronize the operation of different 
components [23]. 

There are two reasons why the frequency of beating 
and self-pulsations do not coincide precisely. Firstly, the 
structure considered in this paper belongs to the open 
systems, and thus the resonant frequencies are not real 
quantities. Depending on the value of the quality fac- 
tor, this can limit accuracy to about 1%. Secondly, the 
Kerr nonlinearity changes the refractive indices of the 
layers and consequently shifts all resonant frequencies. 
This can increase the discrepancy even further, but since 
the period of the self-pulsations is determined by the fre- 
quency splitting of the doublet rather than its absolute 
position, the accuracy of about 10% can be considered as 
acceptable to prove the concept. 

The quantitative agreement can be improved if self- 
pulsations are not very intense. As an example we can 
apply our interpretation of self-pulsations to a paper of 
other authors [H], where this phenomenon was studied 
for microcavities embedded in photonic crystal waveg- 
uides. We derived an explicit formula for the frequency 
of beating in their case and found that it differs from the 
frequency of self-pulsations only in 4%. It is worth not- 
ing that the concept of gap solitons O |4j which is often 



used as a physical explanation of self-pulsations does not 
give any estimations at all in this case. 



V. CONCLUSIONS 

In conclusion, self-pulsations can be viewed as a bal- 
anced state formed by switching and beating of modes 
in coupled nonlinear microcavities. We proposed the de- 
sign of multilayered structure where self-pulsations can 
be excited very efficiently. Since the phenomenological 
CMT equations were applied to describe the switching 
dynamics, the specific choice of configuration is not crit- 
ical, and we expect qualitatively similar results in other 
configurations and physical systems. 

It is also possible to observe self-pulsations in systems 
with continuous spectrum such as waveguide couplers or 
dual-core fibers [551 HHj • The coupling between waveg- 
uides leads to the beating of modes in space, and the 
overall dynamics can be described by equations which 
are similar to Eqs. ([6])-([7]). However, it is not easy to ex- 
cite selectively odd or even modes in that case [27], and 
therefore the usage of coupled microcavities offers some 
advantage. 
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